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ABSTRACT 



PSR J1903+0327 is a mih-second pulsar with a mass of 1.67Mq in a highly eccentric 
orbit (e = 0.44) around a main-sequence star. This unique system cannot be reconciled 
with current observations where milli-second pulsars are generally seen to orbit white 
dwarfs in almost exactly circular orbits. Current theoretical models of binary and stellar 
formation and evolution cannot explain the high eccentricity of this system either. In 
this work, we present three new epochs of optical spectroscopy for the companion to 
PSR J1903-I-0327, obtained to confirm the association of the main-sequence star with the 



pulsar. These 3 new epochs together with the 2 previous ones obtained by lFreire et al. 
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(|201ll ) firmly establish the high eccentricity of the companion's orbit as predicted by 
pulsar timing. Using all five epochs of optical data, we have provided an independent 
estimate of the mass-ratio, R = 1.56 it 0.15 as well as the systemic radial velocity of the 
binary, 7 = 42.1 zb 2.5 km s~^. We constrain the spectral type of the pulsar companion 
to lie between F5 V - GO V (a slightly earlier type than suggested previously) by 
measuring the equivalent widths of two of the three Ca-triplet lines (8498 A, 8542 A) 
and the Paschen line at 8598 A (P14); we also broadly constrain the metallicity of the 
companion. Additionally, we have placed a somewhat better limit on the rotational 
velocity of the pulsar companion of Uj,otsini < 66 km which is still not sufficient to 
allow a test of general relativity using this system. 

Subject headings: stars: neutron, pulsars, stars: pulsars: individual: (PSR J1903+0327) 



1. Introduction 



Binary pulsars are excellent systems for understanding the nature of binary star evolution 
and they also act as laboratories for tests of extreme physics. The milli-second p ulsar (MSP) 



J190 3+0327, was first discove red in the Arecibo L-b and Feed Array pulsar survey (jCordes et al 
20061 ). It was soon found bv IChampion et al.l (|2008l ) that the MSP has a spin period of 2.5 ms 
and is in a highly eccentric orbit (e = 0.44) around what appeared to be a solar- mass companion. 
It is atypical of MSPs to exhibit such high eccentricity or to have a main-sequence (MS) co mpan- 
ion. Current models of stellar and binary evolution do not predict such a scenario either ([Stairs 
20041 ). This is the first and only MSP in the Galactic disk that exhi bits such orbital attrib utes, 
and therein lies its uniqueness. To search for the pulsar companion, IChampion et al.l (120081 ) ob- 
tained near- infrared J-, H- and images of the pulsar field with the Gemini-North telescope and 
found a single star at the position of the pulsar with J=19.22, H=18.41 and Ks=18.03 magnitudes. 
They determined, from the density of stars in the field, that the probability of finding a star at 
the position of the pulsar is < 3%. However, it was not known whether this star was a chance 
encounter or a distant third mein ber in a hierarchical triple system. To confirm this association 
or lack thereof. iFreire et al.l (|201ll ): F2011 hereafter, obtained optical spectroscopy of the possible 
pulsar counterpart using the FOcal Reducer and low dispersion Spectrograph (F0RS2) on the Eu- 
ropean Southern Observatory's (ESO) Very Large Telescope (VLT). They determined the spectral 
type of the companion to lie between an early to mid-G dwarf and from radial velocity changes, 
they positively confirmed the MS star's close-in association with the pulsar. It is typically seen that 
MSPs with spin-periods < 10 ms have whi te dwarfs (WD ) as their companions and exhibit very 
circular orbit s with eccentricit i es, e < .001 (jPhinnev Ill992l ). To explain this inconsistency in PSR 
J1903+0327, IChampion et al.l (j2008l ) proposed the possibility of a triple system where the pulsar 
is orbiting an unseen massive white dwarf (WD) of 0.9-1.1 Mq while the MS star is in a much 
wider orbit around the MSP- WD binary and provides the eccentricity to the inner binary via the 
Kozai mechanism, which causes periodic exchange between inclination and eccentricity of the orbit 
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(|Kozailll962l ). Recently, I Z wart et al.l (|201ll ) have questioned the plausibihty of the existence of such 
a MSP-WD-MS triple system and have pointed out several shortcomings associated with such a 
proposition. It should be mentioned that triple systems with eccentric orbits have been observed 
in the past. For example, the soft X-ray transient, 4U 2129+47 (V1727 Cyg) is a low mass X-ray 
binary (LMXB) with an or bital period of 5.24 hr and a third companion (F-type M S star) in an 



eccentric orbit of 175 days (iGarcia et al 



198S 



Bothwell et al. 



200 



Lin et al.ll2009l ): similarly the 



cataclysmic variable EC 19314-5 915 has an orbital period of 4.85 hr w ith a G8 V star as the outer 
companion (jBucklev et al.lll992l ). An alternative scenario proposed bv lChampion et al.l (|2008l ) was 
the possibility of an exchange interaction in a globular cluster that eventually resulted in the MSP 
orbiting the MS star as its present companion, but the new observations of F2011 have shown this to 
be highly unlikely. Other formation and evolution mechanisms discussed in F2011 are summarized 
in Section 4.2. In the same section, we also discuss the results obtained from the numerical simu- 



lations bv lZwart et al.l (|201lh to test the formation scenario of PSR J1903+0327 proposed by F2011. 



Based on the observations of lChampion et al.l (|2008l ) and F2011, we appear to have the unique 
case of a MSP in orbit around a MS star in a highly eccentric orbit in the Galactic disk. This can 
prove to be of immense importance in the field of general relativity. Pulsar timing of this unique 
system can test alternative theories of gravity and probe the equation of state (EOS) of super- 
dense matter. Accurate pu lsar masses already set interesting limits for the EOS of neutron star 
(NS) material ()Freirel l2009l ) . ruling out a few softer models. Also, the rate of orbital decay of this 
system will constrain various scalar-tensor theories of gravity because the members of this binary 
system have vastly different self-gravities, ( 0.2 for NSs and 10"^''' for MS stars). An accurate 
measurement of the orbital decay of the pulsar companion accompanied by the confirmation that 
it is a slowly rotating solar mass MS star can constrain a significant portion of parameter space 
for general relativity versus scalar tensor theory. So far, the best candidate for constraining these 
dipole contributions consists of a 394 milli-secon d pulsar orbiting a WD with comparable mass 
in a moderately eccentric orbit (jBhat et al.ll2008l ). In F2011, the rotational velocity of the pulsar 
companion could not be very tightly constrained ( < 140 km s~^, 3-0"). An accurate estimate of 
the rotational velocity of the pulsar companion is important to constrain the classical contributions 
to the measured periastron advance of the MSP. The dominant classical contributions arise from 
the rotationally-induced quadrupole moment of the star and this classical spin-orbit c oupling ca n 
account for 10% of the measured periastron advance for a rotational period of a few days (IWexlll998l ) . 



In this paper, we present three new epochs of optical spectroscopy of the pulsar companion 
obtained with the Gemini Multi-Object Spectrograph (GMOS) at Gemini-N. In Section 2 we de- 
scribe the observations and data reduction procedures. In Section 3, we present a detailed analysis 
to confirm the results of F2011 as to the association of the pulsar companion with the MS star, 
determine the spectral type of the pulsar companion as accurately as possible using a combination 
of F2011 data along with our new spectra and constrain the rotational velocity of the pulsar com- 
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panion more precisely. Section 4 discusses the various formation and evolution scenarios possible 
for this system. The result obtained on the rotational speed of the pulsar companion has several 
implications in the field of general relativity which are also presented in this section. Finally, we 
summarize the conclusions in Section 5. 



Observations and data reduction 



Long slit (1.0" slit-width), moderate resolution (R=2200) spectroscopy of the pulsar compan- 
ion m the 7500-9600 A range was performed with GMOS on Gemini North on three different 
observing epochs separated by 20 - 40 days. The observations took place on Aug 19, Sep 06 and 
Oct 05, 2010 using the EEV detector with five 30 minute individual exposures (total of 2.5 hours 
on-source integration time). Wavelength calibration exposures were obtained before and after each 
target exposure, which were averaged to determine a wavelength solution for that particular ex- 
posure. Wavelength solution accuracy was further checked using the sky lines in the spectrum. 
Additionally, a bright slowly rotating (11 kms~^) F-type star wa s observed before and after the 
target observations which served as our rotational standard star (jWolff &: SimonI 119771 ). The de- 
tector was read out with no binning, giving a resolution of 3.4 A sampled at 0.34 A pix~^. The 
Gemini IRAF package for GMOS data reduction was used to reduce the data which involved the 
processes of bias subtraction, flat-fielding, sky-subtraction and finally, spectral extraction. No flux 
calibration was performed on the data. F2011 points out that a bright star located 2.3" from the 
pulsar counterpart contributes partly to the counts from the latter, complicating its ext r action 
in their spectrum. Hence they used the optimal spectral extraction technique by iHvnesI (|2002l ) 
to obtain the spectrum of the pulsar companion. In Fig. 1, we show the position of the target 
compared to the bright source located 2.3" away. Upon fitting Moffat profiles to both the target 
and the bright object, we notice that in the better seeing conditions of the GMOS observations, 
the contamination caused by the bright object to the target is negligible. Additionally, we selected 
a narrow spatial window (0.5" - 0.6") for spectral extraction of the target, thus further reducing 
any contamination from the bright source. Our final extracted spectra still has residuals from 
night sky lines. Furthermore, the wavelength range from 8600 — 8700 A is affected by the broad 
feature in the airglow spectr um which is contributed by a blend of R- and P- branches of 02(0 — 1) 



(Broadfoot Kendall 



19681 ). This, along with residual night sky lines, resulted in our inability to 



recover the Ca II 8662 A feature in the spectrum of the pulsar companion. In the remainder of the 
paper, we will focus on the wavelength region 8400 - 8600 A where most of the prominent features 
reside and w here our sky subtractio n was best. The features in this region were i dentified using the 



line li st from ICenarro et al.l (l2001al ) and the night sky lines were identified from lOsterbrock et al 
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3. Analysis 

In Section 3.1, we confirm the results of F2011 as to the association of tlie MS star to the pulsar 
by comparing the radial velocities obtained by us in three different epochs against the radial velocity 
of the pulsar companion obtained by F2011. Second, in Section 3.2, we constrain the spectral type 
of the pulsar companion by employing equivalent widths of some of the prominent absorption line 
features in our spectrum and comparing them against values from the literature. In Section 3.3, 
we have constrained the rotational velocity of the pulsar companion using the rotational standard 
star that was observed for this purpose. 



3.1. Radial Velocity of the pulsar companion 

To compute the radial velocities at each of the three epochs, we first obtained the wavelength 
shift in the strongest Ca II triplet feature at 8542 A. The accuracy of the wavelength calibration was 
established by comparin g the position of n i ght sk y lines in the spectrum of the pulsar companion 



against the values from lOsterbrock et al.l (jl996l ). From this comparison, we found zero point 



offsets in the wavelength calibration of -0.31 A, -0.25 A and -1-0.15 A for the three observation 
epochs respectively, which were then taken in account while computing the final radial velocities. 
After shifting the individual spectra by their respective wavelength offsets, we created an averaged 
spectrum of the pulsar companion shown in Fig. 2. The calcium 8542 A and the 8498 A features 
are seen in the spectrum. The absorption feature at 8598 A is at the location of the P14 line 
in the Paschen series. These features are present in the spectrum at each epoch. There are also 
possible detections of Fe I at 8514.1 A and 8688 A. However, considering the modest S/N (~ 7 - 
10) of the spectrum, we have decided to restrict our analysis to the lines : Ca II (A8498, A8542) 
and P14(A8598). The deviation of the data points from a straight-line fit to various regions of 
continuum was used to assign the error on the normalized spectrum. In Fig. 2, the location of 
sky line residuals have been replaced with similar continuum regions from the CMOS spectrum 
of the rotational standard. These locations are marked with vertical dashed lines in Fig 2. To 
refine our estimates for the radial velocities, we cross-correlated the averaged spectrum with the 
individual unshifted ones at each epoch while correcting for the wavelength calibration accuracy. 
The cross-correlation function in each case was fitted with a Gaussian and the error was obtained 
by the uncertainty in the determination of the line centroid. This error was then propagated in 
the estimation of the radial velocites. As a result, we obtained radial velocities of Vi = 39.1 it 4.5 
km s~^, V2 = 25.2 ± 5.4 km and V3 = 17.8 zt 6.4 km with respect to the solar system 
barycenter at the three different epochs, respectively. Using the radial velocities obtained by us 
and F2011 at five different epochs, we can estimate the systemic radial velocity, 7 of the binary, 
which cannot be derived from radio timing. Towards that end, we minimized the chi-square, xt^ 
between the observed radial velocities and the corresponding predicted values using the parameters 
determined from radio timing, while varying 7. By using all five data points, we have refined the 
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estimate for the systemic velocity given in F2011 to 7 = 42.1 ± 2.5 km s"^ {xi^ = 0.65, la), a value 
which supports the suggestion of F2011 that PSR J1903+0327 did not originate in a dense stellar 
cluster. Using this value of the systemic velocity, we next minimized the chi-square as before while 
varying the mass-ratio alone. From this, we obtained an independent estimate of the mass-ratio, 
R, and refined the value given in F2011 to i? = 1.56 it 0.15 {xu = 0.50, la). This estimate is in 
agreement with the result derived from radio timing. We did a similar analysis to constrain the 
eccentricity, e, of the binary using five epochs of optical observations and found this number to lie 
between 0.42 < e < 0.58 (xv = 0.40, lo"), thereby confirming the high eccentricity of the binary. 
Further optical/infrared observations of the pulsar companion could be useful in placing tighter 
constraints on the parameters discussed above. Fig. 3 shows the predicted radial velocity curve of 
the pulsar companion (solid line), based on the orbital parameters determined from pulsar timing 
and the systemic radial velocity obtained from five epochs of optical data. The dotted-dashed lines 
in the same figure shows the radial velocity curve of the pulsar companion for an eccentricity of 
e = 0.50, mass-ratio of R = 1.62 and systemic velocity of 7 = 42.1 km s~^ (dotted-dashed line). 
Our observations confirm the result of F2011 that the optical counterpart is indeed the close-in 
companion to the pulsar, orbiting it in a 95-day period. With five epochs of radial velocities 
in-hand and in agreement with the predicted reflex motion of the pulsar companion, the highly 
elliptical orbit of this system is confirmed as well. 



3.2. Spectral type of the pulsar companion 



Significant work has been done in the past to use equivalent widths (EWs) of the calcium triplet 
as well as the Paschen lines as indicators of temperature, metallicity, and luminosity of a star. This 
is the method we have employed in det ermining the s pectral type of the p ul sar companion. T o 
this end, we have follo^ yed the works of MallikI (119941 ). iGinestet et al.l ()l994l ). Ijones et al.l (j 19841 ). 
Idiart et al.l (j 199 71 ) and lDiaz et al.l (jl989l ). The spectral features along with the line windows used 
for calculation of EWs and our results are shown in Tables 1 and 2. EWs were calculated from the 
the averaged, normalized spectrum and errors on EW values were obtained by choosing different 
continuum points for normalization (it lex). Since the Ca II line at 8662 A was not recovered in our 
spectrum, we have used the EW of this feature from the normalized spectrum of F2011 wherever 
applicable (see Table 2). Furthermore, we have reanalyzed the data of F2011 using the two well- 
detected Ca II features at 8542 A and 8662 A shown in Table 3. This was done to ascertain if 
tighter constraints could be placed on the spectral type of the pulsar companion from their data. 
It should be noted that the error in the EWs obtained from their already normalized spectrum is 
based on the error in the Gaussian fit to the feature rather than on continuum placement. 

A summary of our conclusions from various references are shown in Table 4. The sum of EWs 
of the calcium triplet is a robust indication of the luminosity class and following all the references 
cited above, the pulsar companion clearly belongs to a luminosity class 'V. It is a determination 
of the temperature that poses a greater challenge. Using the EWs of the calcium triplet alone does 
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not constra in the spectral type a ny tighter than between F5 - G5 ( see Table 4). On the other 



hand, from iGinestet et al.l (|1994| ). the non zero EW of the P14 line restricts the spectral type to 
GO or earlier. Our confidence in the P14 detection is supported by the fact that we observe this 
feature in all three epochs of data with the same wavelength shift as the calcium lines at 8498 A 
and 8542 A. Taking this into account, we conclude that the spectral type of the pulsar companion 
lies between F5 V - GO V. The last column in Table 4 summarizes the metallicity obtained for 
the pulsar companion from each of the references. While indicative, we do not attempt to draw 
any significant conclusion from it since the metallicity results are limited by the uncertainty in the 
spectral type. 

F2011 determined the spectral type of the pulsar companion by comparing their averaged 
spectrum with synthetic spectra and obtained a temperature of Tgjf = 5825 it 200i^ and a surface 
gravity of > 4. cm s~^. Their analysis points to a star roughly of spectral type GO - G5 V 



( Ali et al.lll995l ). However, from an equivalent width analysis of the Ca II features at 8542 A and 
8662 A obtained from their data, we can conclude only that the star has a spectral type between 
F5 V - G5 V. But we also note that they appeared to have detected P14 as we do although at the 
low signal-to-noise based on our own measurement (see Table 4), constraining the spectral type to 
GO or earlier. There is a slight mismatch between our conclusion and F2011 regarding the spectral 
type of the pulsar companion ( F5 V - GO V vs GO V - G5 V) . However, the detection of Paschen 
P14 line in all three epochs of our data and that of F2011 requires a spectral type GO or earlier. 
Taking both studies into consideration, GO V appears to best describe the spectral type of the 
pulsar companion. This also matches the prediction of the companion's mass of 1.667M0 from 
pulsar timing (F2011). 



3.3. Rotational velocity of the companion star 



We observed a bright rotational standard star (HD 176095) with GMOS at each epoch before 
and after our target obs ervations. The rotat ional standard star is an isolated, slowly rotating (11 
kms~^) late F type star ( Wolff &: Simonlll977l ). Fig. 4 shows the spectrum of the pulsar counterpart 
plotted over the spectrum of the rotational standard (in red). At a spectral resolution of 136 km 
s~^, the effect of broadening from a star rotating at 11 km s^^ cannot be detected. Further, the 
calcium line s in solar-like stars ar e often broadened by collisional or pressure effects in the star's 
atmosphere (jSmith &: Drakd 119871 ). From the spectrum of the pulsar companion, we find that the 
FWHM of the Fe I feature at 8515 A, which is insensitive to pressure broadening, is similar to the 
FWHM of an arc lamp line, indicating that the pulsar companion is not rotating at a speed higher 
than the instrumental resolution. Therefore, an upper limit on the rotational velocity of the pulsar 
companion can be set by finding the minimum broadeni ng that is detectable a t the resolution of 
our in strument using cross-correlation techniques. From iTonrv fc Davis I (119791 ) and iBailer-Jones 
(j2004l ). the measured width (FWHM) of the cross-correlation function (CCF) can be expressed as 
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(^rot + ^nnt + 2o",^ 



inst 



(1) 



where arot is the rotational broadening, anat is the intrinsic broaden ing and amst is the instrumental 
broadening. Following the argument made in Isailer- Jone j ( 2004 ). arot can only be detected if it 
exceeds \/{'2')o'inst, assuming that anat is negligible. Also, the full width of a rotational profile is 
contributed equally by both the approaching and the receding limb of the star and hence corresponds 
to twice the rotational velocity. The oretically, we shou ld be able to detect a minimum rotational 
velocity of Vrotsini = ^J\2){ainst/'^) ( Bailer- Jonesll2004l ). 



To test that, we obtained the autocorrelation function (ACF) of the rotational standard star 
separately in the wavelength ranges 8000 - 8450 A and 8670 - 8800 A . These regions mostly 
comprise the Fe I, Mg I and Ti I features while excluding the pressure sensitive calcium features. 
The average of the FWHMs of these two wavelength regions on either side of the calcium triplet 
was used as a measure of ameas- We then artificially broadened the spectrum of the rotational 
standard in steps of 3 km s~^ and cross-correlated the broadened spectrum with the original un- 
broadened one until a measurable difference was noticed in the CCF. From Eq. 1, cr^gas — "^^Inst 
(for the ACF) and Cmeas = '^'rot + '^'^fnst i^^^ CCF) and subtraction of one from the other gives 
f^rot- We determined FWHMs of the ACF and the CCFs by fitting them with Gaussian profiles. 
A measurable difference in the CCF was obtained when the rotational standard star spectrum was 
broadened by arot = 132 km s~^. Hence, the minimum detectable Vrotsimi at our instrumental 
resolution is 66 km s~^. This value is close to the theoretical minimum de t ectab le rotational speed 
of 96 km s~^ (=136 x -v/2 /2), obtained from the argument in iBailer- Jones (2004). Therefore, these 
observations constrain the pulsar companion's rotational speed to < 66 km s~^. Our constraint 
on the rotational velocity of the pulsar companion is a factor of ~ 2 less than the value quoted in 
F2011. 



4. Discussion 

From our analysis, we have confirmed that the the observed star in the vicinity of the MSP is 
indeed its companion orbiting it in a highly elliptical 95-day orbit as predicted from pulsar timing. 
From the EWs of the calcium triplet and the P14 feature, we constrain the spectral type of the 
star to lie between F5 - GO V. Additionally, we have proven that the companion is rotating slowly 
with an upper limit of <66 km s~^. In this section, we will discuss the implication of our result 
in two areas : constraining evolution scenarios of this unique system and using this system as an 
astrophysical laboratory for tests of general relativity. 
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4.1. Implications for General Relativity 

The present analysis should allow us to constrain the classical spin-orbit contribution to the 
measured rate for the periastron advance: tOo = 86.38 it 0.08 arcsecond century"^. This is because 
the above contribution arises from the rotationally induced quadrupolar moment of the main- 
sequence companion in our b i nary p ulsar and in what follows we provide a few estimates for uJso- 



Using equation (79) in IWexl (jl998l ). we obtain a maximum value for the apsidal motion due to 
classical spin-orbit coupling to be ujso ~ 7.6 x 10^ x J2 arcsecond century"^, where J2 is the 
quadrupolar moment of companion star that we observed. If we let J2 take a value close to that for 
the Sun ( J20 ~ 1-7 x 10"''), we get uso ~ 0.013 arcsecond century"^ and this is roughly an order 
of magnitude smaller than the measurement uncertainty in a;o(F2011). However, the fact that J2 
is proportional to the square of the angular velocity of the star's proper rotation implies that we 
may use J2 ~ </2 x ^^rot/^O' where Vrot stands for the rotational velocity of the companion and we 
let ^0 ~ 2 kms~^. Therefore, we can provide few constraints for ujso contributions to ujo- 

If we let the spin axis of the star be parallel to the orbital angular momentum of the binary, 
we have Vrot ~ 66 km s~^ which leads to ujso ~ 14 arcsecond century"^. This is even substantially 
higher than the 2.3 arcsecond century"^ upper bound on the spin-orbit contributions to the apsidal 
motion that was obtained by treating general relativity to be the correct theory and using the 
minimum total mass compatible with the measurements of three relativistic observables (F2011). 
Interestingly, the fact that we have restricted the spectral type o f the companio n between an F5V 



- GO V, implies that Vrot may be in the range 12 - 25 km (jTassoull l2007l ) . These rotational 
velocity bounds lead to two estimates for ujso- 0.5 and 2 arcsecond century"^ respectively. These 
are also higher than the present measurement error in cOq indicating that it will be rather difficult to 
perform a test for general relativity using the available measurements of the three post-Keplerian 
parameters (F2011). 

Our 8 meter spectra were taken under good to excellent conditions and yet have a relatively 
modest signal to noise at a resolution of R=2200 which is still ~ 7-8 times poorer than what we need 
to achieve a constraint of less than 10 km on Vrot sin i. The lines which are unaffected by pressure 
broadening in the star's atmosphere are too weak to be detected clearly at the current resolution 
and S/N and would require much better signal-to-noise to be confidently detected. Therefore, 
employing the agreement of the three different post-Keplerian parameters to perform a test of 
general relativity is not plausible at the current time. 



4.2. Constraining evolution scenarios for the pulsar system 



The detection of a MS star as the companion to the pulsar and the confirmation of a highly 
elliptical orbit raises vari ous questions on the o rigin of this system while eliminating some other 
hypotheses. For example, IChampion et al.l (j2008l ) considered the possibility of a triple system where 
the pulsar is orbiting an unseen white dwarf (0.9 - I.IMq) in a 95-day orbital period while the 



- 10 - 



MS star ha ppens to be in a longer period that drives the eccentricity of the system by the Kozai 
mechanism (|Kozailll962l ). Since we have estabhshed that the MS star orbi ts the pulsar in a 9 5-day 
orbital period, the triple system hypothesis can be ruled out. Additionally, I Z wart et alj (|201ll ) have 
pointed out that the hypothesis of a MSP-WD-MS triple system is precluded by the fact that the 
observed change in the eccentricity of this system is three o rders of magnitude smaller than that 
predicted via the Kozai mechanism (jGopakumar et al.ll2009l ) . 



Champion et al.1 (|2008l ) also discuss the possibility of an exchange interaction in a dense stellar 



environment that might have brought the MS star to orbit the pulsar in a 95-day period. In F2011, 
the probability of this was found to be negligible. One of the parameters used in their simulation is 
the systemic radial velocity, 7, of the binary. With five epochs of radial velocity data, we have now 
refined this number to be 7 = 42.1±2.5 km s~^. This is consistent with the value obtained by F2011 
thereby supporting their suggestion that the system did not form in a dense stellar environment. 

F2011 also find that the pulsar has been recycled or spun- up to milli-second speeds (due to 
mass-transfer from a past, accreting X-ray binary phase) but show that the current MS star is 
unlikely to have played the role of the donor star responsible for its recycling. To answer the 
question of how the pulsar was spun-up and what happened to the donor star, F2011 investigated 
another triple system scenario. They considered a triple system wherein the outer companion (the 
present MS star) started out in a much wider orbit around the inner binary, which consisted of 
two more massive MS stars in a shorter orbital period. The more massive star in the inner binary 
evolves into a red supergiant phase and then engulfs the companion into a common envelope phase. 
As the orbit of the inner binary decays, it will release orbital energy sufficient enough to leave a 
much closer binary consisting of a He-core and an almost unaffected MS star. Meanwhile, if the 
outer star in the triple system is close enough, it could be engulfed by the expanding envelope of 
the inner binary as well, which will cause it to spiral in and result in a much closer triple system. 
Once the He-core explodes in a supernova, we are left with an inner Low-Mass/Intermediate-Mass 
X-ray Binary (LMXB/IMXB) and an outer MS star. The NS is then spun up to milli-second 
periods by the donor star, which has filled its Roche-lobe by that time. As mass transfer continues 
from the donor to the NS, the orbit of the inner binary in the triple will widen to some extent 
and will end in a chaotic three-body interaction. This will result in expulsion of the less massive 
star from the system, eventually leaving the MSP in a tighter and eccentric orbit around the MS 
star (the initial outer body in the triple system). F2011 have also discussed in less detail several 
alternative scenarios for the formation and evolution of the system. All of these scenarios start 
out with the system as part of a triple and then invoke different dynamics that transpires after 
the common envelope phase, leading to the system as we see it now. However, thorough numerical 
simulations are necessary to provide more insight into the formation and evolution of these systems 
in a non-deterministic way. 



Recently, IZwart et al.l (j201ll ) performed a detailed numerical analysis to test the formation 
scenario presented in F2011. They start out with a stable LMXB (Porb < l.Od) with a MS star in 
a wide and possibly eccentric outer orbit. They have shown from their simulations that as mass 
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transfer starts in the inner binary and as the orbit of the inner binary widens, the system becomes 
unstable. Depending on the exact configuration of the system at this point, each one of the three 
components in this triple system could be ejected. Specifically, they found that a stable hierarchical 
triple consisting of an inner binary with component masses of 9-13 Mq and 0.8 - 2.0 Mq, separated 
by > 2OOi?0 and an outer MS star with mass < 2Mq and having a semi-major axis > 560i?Q, could 
eventually result in the formation of a system like PSR J1903-I-0327. Even though the chances of 
passing through the vario us chain of events w ere found to be low, the result matches with what we 
see in PSR J1903-I-0327. I Z wart et al.l (|201ll ) were also able to describe an evolution scenario for 
the soft X-ray transient 4U 2129+47 (V1727 Cyg) and the cataclysmic variable, EC 19314-5915, we 
noted before in Section 1. Additionally, they concluded that systems like PSR J 1903+0327 have a 
birthrate less than 3% compared to the Galactic LMXB's and 10 times less than all MSPs. Thus 
the probability for this formation scenario happening, is not large, but it is also not nil. Their 
model predicts tens to hundreds of systems similar to PSR J1903+0327 existing in our Galaxy and 
raises the question of why no system like this has been found prior to PSR J1903+0327. 



In another independent work. iLiu fc Lil (120091 ) consider that PSR J1903+0327 started out as a 
binary system but during the supernova (SN) explosion that produced the pulsar, a fall-back disk 
was formed around the newborn NS. Accretion started from the disk onto the NS, which spun it 
up to milli-second speeds. Meanwhil e the high eccen tricity is considered to be a consequence of the 
SN explosion that produced the NS. iLiu &: Lil ((20091) mention that if this system was undisrupted. 



it could maintain this eccentricity for > 10 yr. One of the uncertainties mentioned in their 
work is whether the fall-back material would have enough angular momentum to produce the disk. 
This calls for an unusual SN-fallback history to explain the very short period (2.5 ms) of PSR 
J1903+0327. In order for the fallback disk scenario to work, it requires the companion of the 
pulsar to be < (1 — 2) X 10^ yr. From our work as well as from F2011, we find the MS star to be 
older than what is required for the fallback disk scenario to work. 



5. Conclusion 



By obtaining new spectroscopy with GMOS at Gemini-N, we have confirmed the results of 
F2011 which positively identified a MS star as the close-in pulsar companion. Radial velocities are 
now in-hand for 5 epochs and the predicted radial velocity curve obtained from pulsar timing is 
confirmed. The highly eccentric orbit is confirmed as well. Additionally, we obtained estimates for 
the mass-ratio of i? = 1.56 ± 0.15 and the systemic radial velocity of the binary of 7 = 42.1 ± 2.5 
km s~^. These refined estimates are important in the context of explaining the origin of this 
system. We have determined the spectral type of the pulsar companion to lie quite close to GO 
V, from an analysis with our data and a reanalysis of F2011 data. Also, we have constrained 
the rotational velocity of the pulsar companion as < 66 kms~^. With the present estimate of the 
Vrotsini, it will be difficult to perform a test for general relativity. Radio pulsar timing can lead 
to better constraints on General Relativity in comparison to scalar tensor theories and should be 
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pursued over the next few years. We have also noted that the evolutionary scenario possible for 
the formation and evolution of th is unique pulsar system is best described by the recent numerical 
simulations of I Z wart et alj (|201ll ). Starting with a stable hierarchical triple system consisting of 
an inner binary with component masses of 9-13 Mq and 0.8 - 2.0 Mq, separated by > 2OOi?0 and 
an outer main sequence star with mass < 2Mq and having a semi-major axis > 56O-R0, numerical 
simulations were able to replicate the formation and evolution of a system like PSR J1903+0327. 
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Table 1. Equivalent Width Analysis 



Feature 


Integration 


limits (A) 


Equivalent width(A) 


CaT 8498 


8495 - 


- 8501 


0.61 ± 0.27 


CaT 8542 


8535 - 


- 8550 


1.70 ± 0.74 



Note . — Tabic of EWs for comparison with the work by 
iMallikI l ll994h 



Table 2. Equivalent Width Analysis 



Feature 


Integration 


limits (A) 


Equivalent width(A) 


CaT 8498 


8488 - 


- 8510 


1.17 ± 0.50 


CaT 8542 


8530 - 


- 8556 


2.32 ± 0.70 


CaT 8662 


8652 - 


- 8672 


1.36 ± 0.15^^ 


P14 8598 


8585 - 


- 8605 


.70 ± .50 



Note. — T a ble o f EWs for compariso n with the work by 
iGinestet et al.l lll994l'l. Ijones et al.l lll984l) . Ildiart et all |l993) 
and iGarcia- Vargas et al .1 l|l998h 

''Da ta for this absorption feature was used from iFreire et al.l 
1I2OIII) 



Table 3. Equivalent Width Analysis 



Feature 


Integration limits (A) 


Equivalent width(A) 


CaT 8542 


8527 - 8557 


2.40 ± 0.17 


CaT 8662 


8647 - 8677 


1.32 ± 0.15 


P14 8598 


8585 - 8605 


.11 ± .07 



Note. — Tab le of EWs for comparison with the work by 
iDiaz et al. I lll989h . The line window for P14 is used from Table 
2. All other line windows arc from lDiaz et al.l \l98^ . The data 
for computa t ion o f the above EWs is obtained from Figure 2 of 
iFreire et al ] 1I2OI1I ') 
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Table 4. Summary of Results 



Reference 


Feature used 


Luminosity Glass 


[Fe/H] 


Spectral Type 


Mallik (1994) 


CaT{A8498+A8542) 


IV/V 


-0.3 < [Fe/H] <0.2 


F5 - G5 


Ginestet et al. ("1994") 


CaII(A8542), CaII(A8542), P14 (A8498) 






F2 - GO 


Ginestet et al. (1994) 


GaT(A8498+A8542+A8662) 


V 




a 


Jones et al. (1984) 


GaT(A8498+A8542+A8662) 


V 




F5 - G?'' 


Idiart et al. ( 1997) 


CaT(A8498+A8542+A8662) 




-0.6< [Fe/H] <0.2 


a 


Diaz et al. (1989") 


CaT(A8542+A8662) 


V 


-0.4 < [Fe/H] <0.2 


F5 - Gs'' 



Note. — Summary of results using the GaT features from our data and F2011 

''A8662 feature was used from F2011 

''both features were used from the data of F2011 
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Fig. 1. — The position of the pulsar (hne # 2442) and the bright object (Une # 2410) 2.3" from the 
pulsar is shown in the above figure for a total of 100 columns. Moffat fits to their profiles indicate 
negligible contamination to the counts from the pulsar companion. 
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Fig. 2. — The normalized averaged spectrum of the pulsar companion is shown with representative 
error bars. Ca II (A8498 and A8542) lines are visible in the spectrum as well as the P14 line at 
8598A and Fe I at 8515A. There are also possible detections of Ti I (8435 A), P17 (8467 A) and 
Fe I (8688.5 A) lines but we cannot be certain due to the low S/N of the spectrum. The position of 
the sky lines are indicated by the vertical dashed lines. Sky lines that are too close in wavelengths 
to be resolved as separate lines, appear as thick dark vertical lines 
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Fig. 3. — The barycentric radial velocity curve of the pulsar (dashed line) and the pulsar companion 
(solid line), as predicted from pulsar timing (e = 0.44, R = 1.62) with 7=42. 1km s~^(obtained from 
optical data). The bottom x-axis and the top x-axis are the Mean Julian Date (MJD) of observations 
conducted by us and F2011 respectively. The data points in filled circles are the measured values 
of radial velocity obtained from our dataset while the two points in filled squares (near MJD 54640 
and MJD 54700) from F2011 are also shown in the figure. The dotted-dashed line represents the 
radial velocity curve obtained for e = 0.50, R = 1.62 and 7 = 42.1 km s~^. 
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Fig. 4. — The normalized and averaged spectrum of the pulsar companion is shown in the figure 
compared to the normalized rotational velocity standard star spectrum (in red). The location of 
the sky lines are marked in vertical dashed lines. 



